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We review the extension of the factorization formahsm for perturbative QCD to 
. soft initial- and final-state scattering associated with hard processes in nuclei. 
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^ : 1 Introduction 

^ ! 

\0 • In this talk, we would like to review a few results from a perturbative QCD (pQCD) treat- 
^ . ment of the scattering of hadrons and leptons in nuclei, based on factorization, work in 
^ ^! collaboration with Ma Luo H, H and, more recently, Xiaofeng Guo 0. At the outset, 
I I it may be useful to clarify the relation of this work to the recent papers of Baier et al. 
(BDMPS), described by Dokshitzer at this workshop [|]. We have tried to illustrate this 
• relation schematically in Fig. |I[ The BDMPS analysis begins (Fig. la) with the classic 
> . treatment of radiation induced when a charged particle passes through a large target, due 
^ , originally to Landau, Pomeranchuk and Migdal (LPM). This analysis does not require the 
^ I presence of a hard scattering, but describes the coherent results of many soft scatterings. 
■ ■ ■ ' Its primary subject has traditionally been induced energy loss. Our analysis (GLQS) begins 
with the perturbative QCD treatment of hard-scattering in a small target (Fig. lb), in which 
the primary subject of interest is momentum transfer. A complete analysis (Fig. Ic) of hard 
scattering in a large target, involves both energy loss and the transverse momenta due to 
initial- and final-state soft scatterings. Our work is a step in this direction, attempting to 
stay as close as possible to the pQCD formalism, in which we may readily quantify cor- 
rections. To be specific, we consider only a single soft inital- or final-state interaction in 
addition to the hard scattering. Our central observation is that for suitably-defined jet and 
related inclusive cross sections this is the first order in an expansion in the quantity 
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where A represents a nonperturbative scale, which we shall identify with a higher-twist parton 
distribution below. That additional scatterings are suppressed by factors of is perhaps 
surprising. Let us review why this is the case, at least for certain cross sections. 




(c) 

Figure 1: Alternate approaches to hard scattering in nuclei, (a) Landau- Pomeranchuk- 
Migdal analysis treats energy loss due to many soft scatterings, (b) Perturbative QCD 
analysis treats momentum transfer due to hard scattering, (c) For scattering in nuclei, both 
must be combined. 

The basic analysis of hard-scattering in nuclear matter (cold or hot) |^ is quite simple. To 
be specific, consider the scattering of a quark. A hard-scattering with momentum transfer 
Q can resolve states whose lifetimes are as short as 1/Q, for instance quarks off-shell by 
order Q, but still less that Q. The off-shellness of the scattered quark increases with the 
momentum transfer simply because the number of available states increases with increasing 
momentum. Similarly, the scattered quark, of momentum p' is typically off-shell by order 
mj < Q. We may think of mj as the momentum of the jet into which quark fragments. If 
we are to recognize the jet, we must have mj <^ Ej = p'^. On the other hand, the counting 
of available states ensures that nij ^ Aqcd- 

Now the scattered quark has a lifetime in its own rest frame At^^ ^ ~ ^ with mj -C 
Ej. In the target rest frame, however, this becomes, for large enough Ej/mj, At^^'^^^^^^ ~ 
m7 (mj) ^ ^^^'^^ (fixed) target size. Thus, at high enough energy the lifetime 

of the scattered quark will exceed the target size, even though the quark itself is far off the 
mass shell, typically by a scale that grows with the momentum transfer Q. 

Further couplings of the off-shell quark are suppressed, first of all by the strong coupling 
evaluated at scale mj, and, more importantly, by an overall factor of 1/m'j ~ since 
the effective size of the scattered quark decreases with momentum transfer in this manner. 

In summary, for inclusive processes such as jet production, high-Q implies that process- 
dependent multiple scattering is power-suppressed compared to single scattering. Initial- 
state interactions internal to the nucleus are leading-power, but factorize. Thus the "Cronin 
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effect", A"-dependence witli a > 1, due to multiple scattering, is liiglier- twist for inclu- 
sive distributions, while the "EMC" effect for parton distributions in nuclei is (almost by 
definition) leading- twist. 

The most important point here is that the scattered particle remains off-shell for its entire 
transit of the target. Thus, its interactions with the target may be treated by the formalism 
of perturbative QCD, which, however, must be extended to include corrections that decrease 
with extra powers of momentum transfer. Up to the first such "higher-twist" contribution, 
a general cross section has the representation 

Cr(Q) = if°®/2®/2+ (^^)if^®/2®/4 + 0(^^^ , (2) 

where ® represents covolutions in fractional momenta carried by partons, and /„ represents 
a parton distribution of twist n. Target-size dependence due to multiple scattering can only 
appear in the second term in this expansion. 



2 Parton-Nucleus Scattering in Perturbative QCD 

2.1 Factorization at Leading and Nonleading Powers 

Let us review some of the details of a factorized cross section like (^. The first term, 
consisting of only twist-two matrix elements has the detailed form. 
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' ^z^i =Y. j dyfj/p2{y,Q) J dxfi/p,{x,Q)aij{xpi,yp2,p') , (3) 

where we may take p' as the momentum of an observed jet. The fragmentation of a jet, 
suitably defined, is calculable in perturbation theory, and may be absorbed into the "hard 
scattering function" a. The fa/p are distributions of parton type a in hadron p. They have 
the interpretaion of expectation values in the hadronic state of products of fields on the light 
cone, for instance, for a quark distribution 

f,/p{x,Q) = J ^e"^"^-(p|g(0)^g(r)b) , (4) 

where for simplicity we choose the ^4+ = gauge, assuming p is in the plus direction. Eq. (j^) 
is ilustrated by Fig. 2a. As shown, the convolution in eq. (^) is in terms of the momentum 
fractions x and y carried by partons i and j, from hadrons pi and p2, respectively, into the 
hard scattering. 

Fig. 2b is the corresponding picture for a higher-twist contribution to hard scattering. In 
this case two partons i and i' with momenta xipi and X2P1 from the target (the "nucleus") 
collide with a single parton j of momentum yp2 (from the "projectile"). 
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Yl J d'yfj/p2iy,Q) j dxidx2dx^T(,i,)ip^{xi,X2,xz,Q)a^^l)^j{xiPuyp2,p') . (5) 
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(a) 



(b) 



Figure 2: Perturbative QCD at leading twist (a), and higher twist (b). 



The expectation value T corresponding to this multiparton contribution from the target is 
typically of the form 0, 



where Bi is the field corresponding to a parton of type i = q,q, G. In eq. (||), the hard part 
depends on the identities and momentum fractions of the incoming partons, but is 
otherwise independent of the structure - in particular the size - of the target (and projectile). 
To find y4-enhancement due to multiple scattering, we must look elsewhere. 

2.2 A- Enhancement from Matrix Elements 

For definiteness, we consider photoproduction or deeply inelastic scattering on a nucleus 
|1|, ^. In this case, the additional soft scattering is always a final-state interaction. The 
structure of the target is manifest only in the matrix element T in eq. (^). Each pair of fields 
in the matrix element (^) represents a parton that participates in the hard scattering. The 
y~ integrals parameterize the distance between the positions of these particles along the path 
of the outgoing scattered quark. In eq. (H), integrals over the distances generally cannot 
grow with the size of the target because of oscillations of the exponential factors e*^^^'^* . 
Poles from a in the Xi integrals, associated with the scattered particle, however, can result 
in finite contributions from points where two of the Xj vanish [1-3]. An example is shown 
in Fig. 1^. It is important to emphasize that using a pole in the complex xt (longitudinal 
momentum) space to do the integral does not correspond to assuming on-shell propagation for 
the scattered quark. Indeed, the Xi integrals are not pinched between coalescing singularities 
at that point, and the same results could be derived by performing the Xi integrals without 
ever going through the Xi = points. 

The result of this reasoning is that matrix elements that depend on three fractional 
momenta, as in (P) above, simplify to a form like 



(6) 
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Figure 3: Pole that gives rise to an A-enhanced cross section. 

X I {vA\mi^F^^{y2)F'-.{y~)q{yr)\PA) , (7) 

where \pa) is the relevant nuclear state. In this form, integrals over the y~ can grow with 
the nuclear radius as fast as A^/^^ once local color confinement is taken into account. The 
variable x here is the fractional momentum associated with the hard parton from the target 
that initiates the process. The soft scattering contributes a negligible longitudinal fractional 
momentum. Details of the reasoning and calculation for deeply inelastic scattering are given 
in Ref. §. 



3 Applications 

In Refs. [|T| and , we have applied the formalism sketched above to single-particle inclusive 
and single-jet production for deeply inelastic scattering and photoproduction. These cases 
involve final-state interactions only. In each case, the leading l/Q^ correction is proportional 
to the matrix element in eq. (0), or to a corresponding matrix element Tq with four gluon 
fields. Of course, the value of the correction cannot be estimated without an idea of the 
magnitudes of the T's. Since these magnitudes are nonperturbative they must be taken 
from experiment. At the same time, we expect the x-dependence of the probability to detect 
the hard parton to be essentially unaffected by the presence or absence of an additional soft 
scattering. Thus, we choose ansatz 

T,{x,Q) = \^A^I^U,,,{x,Q) (8) 

for Tq, in terms of the corresponding twist-two parton distribution /, with A a constant with 
dimensions of mass (see eq. ([I|)). This assumption facilitates the comparison to data. 

A quantity that is sensitive to final-state rescattering in a particularly direct way is the 
momentum imbalance of di-jets in photoproduction in nuclei. The A^^^ dependence of this 
quantity is related to the matrix elements Tq and Tq by the simple formula p| 

(4^^^)4/3= E /dxr,(x,A)i7^"(xp,p„£) = AM^/3 E fdxfaix,A)H'"^ixp,p,,i), 

^ ^ a=q,g-^ a=q,g 

(9) 
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where H"''^ is a hard-scattering function that we have computed to lowest order and where in 
the second equahty we have used (|]). The momentum i may be identified as the momentum 
of the more energetic jet. By comparing eq. (H) to data, [H we found ~ 0.05 — 0.1 GeV^ 
1^. This value may be used to predict anomalous A-enhancement for other processes. 

One such process is direct photon production at measured transverse momentum, whose 
very moderate A-dependence has been measured by the E706 experiment at Fermilab. In 
Ref. 0, it was found that the value of above, which produces a relatively large enhance- 
ment in dijet momentum imbalance, due to final-state interactions, produces a quite small 
^-enhancement in photoproduction, due to initial-state interactions, consistent with exper- 
iment. This may shed some light on the long-standing observation that (initial-state) trans- 
verse momentum effects in Drell-Yan cross sections are also surprisingly small 0. Clearly, 
further study of this and related questions is in order. 
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